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Abstract 

In the present paper the results obtained in the investigation of possible diur¬ 
nal effects for low-energy single-hit scintillation events of DAMA/LIBRA-phasel 
(1.04 ton X yr exposure) have been analysed in terms of an effect expected in 
case of Dark Matter (DM) candidates inducing nuclear recoils and having high 
cross-section with ordinary matter, which implies low DM local density in order 
to fulfill the DAM A/LIBRA DM annual modulation results. This effect is due 
to the different Earth depths crossed by those DM candidates during the sidereal 
day. 


Keywords: Scintillation detectors, elementary particle processes. Dark 
Matter 

PACS numbers: 29.40.Me - Scintillation detectors; 95.30.Cq- Elemen¬ 
tary particle processes; 95.35.-|-d - Dark matter (stellar, interstellar, galac¬ 
tic, and cosmological). 
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1 Introduction 


The present DAMA/LIBRA experiment [1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 
13], as the former DAMA/Nal [14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 
26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43] has the 
main aim to investigate the presence of DM particles in the galactic halo 
by exploiting the model-independent DM annual modulation signature 
(originally suggested in Refs. [44, 45]). In particular, they have cumula¬ 
tively reached a model independent evidence at 9.3(T C.L. for the presence 
of DM particles in the galactic halo by exploiting the DM annual mod¬ 
ulation signature [4]. Recently the results obtained by investigating the 
presence of possible diurnal variation in the low-energy single-hit scintilla¬ 
tion events collected by DAMA/LIBRA-phasel (1.04 ton x yr exposure) 
have been released and analysed in terms of a DM second order model- 
independent effect due to the Earth diurnal rotation around its axis [12]. 
In particular, the data were analysed using the sidereal time referred to 
Greenwich, often called GMST. No diurnal variation with sidereal time 
has been observed at the reached level of sensitivity, which was not yet 
adequate to point out the effect searched for there. In the present paper 
those experimental data are analysed in terms of an effect - named Earth 
Shadow Effect' - which could be expected for DM candidate particles in¬ 
ducing nuclear recoils; this effect would be induced by the variation ~ 
during the day - of the Earth thickness crossed by the DM particle in 
order to reach the experimental set-up. It is worth noting that a similar 
effect can be pointed out only for candidates with high cross-section with 
ordinary matter, which implies low DM local density in order to fulhll 
the DAMA/LIBRA DM annual modulation results. Such DM candidates 
could get trapped in substantial quantities in the Earth’s core; in this 
case they could annihilate and produce secondary particles (e.g. neutri¬ 
nos) and/or they could carry thermal energy away from the core, giving 
potentiality to further investigate them. 

Preliminary investigations on DM candidates inducing diurnal varia¬ 
tion were performed in Refs. [46, 47, 23] and more recently in Ref. [48]. 

2 The Earth Shadow Effect 

During a sidereal day the Earth shields a terrestrial detector with a vary¬ 
ing thickness, and this induces a variation of the flux of the DM candidates 
impinging the detector, mainly because of the modihcation of their veloc- 
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ity distribution, f{v). It is worth noting that this Earth Shadow Effect is 
very small and could be detectable only in case of candidates with high 
cross-section with ordinary matter (i.e. present in the galactic halo with 
small abnndance). 

The detector (and the hosting laboratory) velocity in the Galactic 
frame can be written as: 

Vlabiff Vl^SR T ^0 T ^reti(t) T Vrotifyi (1) 

where: (i) vlsr is the velocity of the Local Standard of Rest (LSR) be¬ 
cause of the rotation of the Galaxy; (ii) Vq is the Snn pecnliar velocity 
with respect to LSR; (iii) Vrevif) is the velocity of the orbital motion of 
the Earth aronnd the Snn and (iv) Vrot{t) is the velocity of the rotation of 
the Earth aronnd its axis. The two latter terms change as fnnction of the 
sidereal time, t. Using the galactic coordinate frame (that is x axis to¬ 
wards the galactic center, y axis following the rotation of the Galaxy and 
the axis towards the galactic North pole), one gets: vlsr = ( 0 , no , 0 ), 
where Vq = (220 ±50) km/s (nncertainty at 90% G.L.) [22, 49] is the local 
velocity, and Vq = (9,12, 7) km/s [50]. 

The DM particles in the galactic halo have a velocity distribntion g{w), 
which depends on the considered galactic halo model. Ref. [22] has shown 
many possible scenarios for the galactic halo; in the following we consider 
the isothermal halo model just because of its simplicity: 

_»2 

g{w) = Ae ~^e{vesc-\w\), (2) 

with A normalization constant and Vesc escape velocity, assnmed in the 
following eqnal to 650 km/s, as often considered in literatnre; however, 
it is also affected by nncertainty. However, no sizeable differences are 
observed in the ontcome when a different valne of Vesc = 550 km/s is 
considered, more closer to the 90% G.L. range of the RAVE Survey results 
[51]. In the laboratory frame the DM velocity distribntion f{v) is obtained 
from eq. 2 straightforward since w = v + viab- 

To evalnate the expected daily variation of the DM particles velocity 
distribntion dne to the Earth Shadow Effect, it is necessary to estimate 
the time dependence of the 9 angle, the “zenith distance” of viab (i-e. the 
distance between viab and the zenith, see Fig. 1). This can be determined 
by astrophysical considerations studying the Earth’s rotation aronnd its 
axis. 

The simplest way to calcnlate 9(t) is in the eqnatorial coordinate sys¬ 
tem [12] where the axis is directed towards the vernal eqninox, and 
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Figure 1: Schematic view of the DM particles impinging direction on a detector; the 
x',y',z' represent the laboratory frame coordinate system. Left: schematic repre¬ 
sentation of the correlation between the thickness, d, crossed by the considered DM 
candidates to reach a laboratory (hypotetically placed at the geographic North pole) 
and the DM impinging angle, Bin- Right: schematic representation of the experimental 
condition considered in the text: detector placed at the Gran Sasso National Labora¬ 
tory (LNGS) with the z' axis in the vertical direction and the x' axis pointing to the 
vernal equinox. 


and 63 “ are on the eqnatorial plane; the 63 “ axis is towards the North 
pole. The right-handed convention is nsed. To work ont the galactic coor¬ 
dinates of those versors, one considers: i) the eqnatorial coordinates of the 
galactic North pole: RA = 192°.859508 and DE = 27°. 128336, where RA 
is the right ascension and DE is the declination; ii) the eqnatorial coordi¬ 
nates of the galactic center: RA = 266°.405100 and DE = —28°.936175, 
evalnated at the Epoch J2000.0. In the galactic coordinates, those versors 
can be written as: 

ef" = (-0.05487,0.49411,-0.86767) 

= (-0.87344,-0.44483,-0.19808) (3) 

= (-0.48384,0.74698,0.45599). 

We dehne Vg = vlsr + %• In this section, when a numerical calculation 
is employed, we assume Vq = 220 km/s; hence, Vg = 232.28 km/s, and in 
the equatorial coordinate system: decs = 42°.18 (the “zenith distance”) 
and (pecs = —46°. 14 (azimuth angle from 6 ®“). For simplicity in the 
following of this section in order to offer an estimate of the Earth Shadow 
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Effect diurnal behaviour we consider Vrev{t) equal to its annual mean 
value, i.e. zero. To introduce the Earth motion around its axis, hrstly we 
consider the simplest case of a laboratory at North Pole and we dehne 
the horizontal coordinate system with z' axis directed as shown in Fig.l, 
left and x' axis directed towards a given longitude Aq. In this coordinate 
system (TV label) the velocity of the Earth can be written as: = Vg, 

On = Oecs, and lpn = <^ecs - + to) = -(v^o + + to)), where 

(po = —‘fees, ‘^rot = /Td with Trf = 1 sidereal day, t sidereal time 
referred to Greenwich, and to = 24Ao/27r sidereal hours. 

The general case of viab in a laboratory at latitude 0o can be derived 
by rotating counterclockwise around y' of an angle a = n/2 — fo: 



= 



( 4 ) 


( COS a 

0 

— sin a \ 

/ 

Vs sin Oecs COs((po + i^rotit + to)) \ 

= 

0 

1 

0 


-Vs sin Oecs sin((^o + (^rotit + to)) 


sino; 

0 

cos a j 

V 

Vs cos Oecs j 


Thus, 0(t) for a laboratory position, identihed by the longitude Aq and 
latitude fo, can be derived from cos 9(t) = ■ z' (see Fig. 1, right) 

obtaining: 

cos 9(t) = sin 9 ecs cos fo COs{Urot(t + to) + fo) + cos Oecs sin 00 ( 5 ) 

= cos 0 cos 00 COs{Urot(t + to) +9:^0) + sin 0 sin 0o, 



GMST (h) 


Figure 2: The angle 0 as a function of the sidereal time in the case of the LNGS 
(latitude Aq = 13°34' E, longitude 0o = 42°27' N). The Earth shielding is maximum 
about at 8:00 h and minimum around 20:00 h; see text. 
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where tjj = 7r/2 — decs = 47°.82 is the altitude. In case of the LNGS, the 
longitude is Aq = 13°34' E and the latitude is 0o = 42°27' N, thus to = 
0.904444 sidereal hours. The 9{t) behaviour at LNGS is shown in Fig 2. 
Note that, before detection, the DM particles at LNGS preferentially cross 
an Earth maximum thickness at about 08:00 h and an Earth minimum 
thickness at about 20:00 h (GMST in both cases). 

3 Deformation of the DM velocity distribution due 

to the Earth Shadow Effect 

To study the experimental data in terms of possible Earth Shadow Effect, 
a Montecarlo code has been developed to simulate the propagation of the 
DM candidates elastically scattering off Earth nuclei in their travel in the 
Earth towards the underground experimental site. For such a purpose 
useful information has been gathered about the Earth composition and 
density. The Montecarlo code numerically estimates the velocity distribu¬ 
tion - in the laboratory coordinate system - of the impinging DM particles 
after having crossed the Earth; such velocity distribution depends on the 
mass of the DM candidate, on its cross-section on nucleons, on the ini¬ 
tial unperturbed velocity distribution, on the sidereal time, and on the 
latitude and longitude of the laboratory: fiab{v,t\mDM,o'n)- Then, this 
velocity distribution has been used to evaluate - in an assumed framework 
- the expected counting rate as a function of the sidereal time in order to 
be compared with the experimental data. 

In this section details are given about the assumptions adopted in 
the simulation, as in particular: the Earth model, the mean free path 
and path reconstruction of such DM candidate, the adopted interaction 
model (nuclear form factor, scaling law, etc.) and the fiab{r,t\mDM,an) 
estimation. 

The Earth model 

An Earth model has been assumed in order to estimate the signal variation 
due to the Earth Shadow Effect] in particular, the matter density and 
composition of the Earth have to be considered. The simulation adopts a 
simplihed Earth model starting by the Preliminary Reference Earth Model 
(PREM) [52]. Just three main Earth layers with a constant density and 
homogeneous distribution in each one (these values are averaged over the 
PREM density distribution behaviour) are considered: i) the Inner Gore; 
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ii) the External Core; iii) the Mantle. The densities and mass percentage 
for each layer are given in Table 1; in this simplihed model the rare 


Table 1: Density values, pL, and t-th nucleus mass percentage, Si, adopted in the 
present calculations for the layers of the considered Earth model [52], 


Layer (L) 

R-mm R-max 

(km) 

Mass percentage 
(Si) 

Density (pi) 
(kg/m^) 

Inner Core 

0 - 1221.5 

Fe (79 %) 

Ni (21 %) 

12839 

External Core 

1221.5 - 3480 

Fe (86 %) 

S (12 %) 

Ni (2 %) 

10901 

Mantle 

3480 - 6371 

0 (44.9 %) 

Si (21.6 %) 

Mg (22.8 %) 

Fe (5.8 %) 

Ca (2.3 %) 

A1 (2.2 %) 

Na (0.4 %) 

4605 


isotopes (mass percentage lower than 0.1%) have been neglected. 


Interactions of the considered DM candidates and path reconstruction 

We assume that the considered DM candidates loss their energy elastically 
scattering off nuclei with spin-independent coupling. The mean free path 
in the L-th Earth’s layer is: 

V = -, (6) 

O'DM,ini 

i=\ 

where: i) Ui = ^5i is the number density of the i-th nuclei in the Earth’s 
layer composed by N nuclear species; ii) <7DM,i is the cross-section on 
the t-th nucleus and rrij is the mass of the i-th nucleus; iii) Pl is the 
layer’s density; iv) Si is the t-th nucleus mass percentage in the layer. We 
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considered a coherent scattering and a scaling law^: 


0'DM,i 


a, 




2/1 


( 7 ) 


where: i) Ai is the mass number of nucleus i] ii) is the DM candidate- 
nucleon cross-section; iii) fj,i is the DM candidate - nucleus(nucleon) 
reduced mass. Thus, assuming nii ~ AiUin, one can write: 


1 

Al 


^nPL ^3 ^ / 'I^DM T 

* * \mDM + AiUln 


2 


( 8 ) 


For simplicity, the deflection of the DM particles crossing the Earth 
is neglected (i.e. a linear trajectory is considered); in this assumption 
the Earth thickness crossed by those DM candidates, d, depends only on 
the impinging angle with respect to the detector, 9in, according to the 
relation: 


d = 2R(^ cos {9in ), (9) 

with Earth radius. This d value is the sum of the distances passed 
through by the DM candidate in each layer: d = J2 l ^l- Dehning the 
maximum radii of the Earth layers Ric, Rec and Rm = R® for Inner Core, 
External Core and Mantle respectively (see Table 1), the number of layers 
crossed by DM particles in the considered schema is: 0 for 9in > 90°; 1 
for 33°.11 < 9in < 90°; 3 for 11°.05 < 9in < 33°.11; 5 for 0° < 9in < 
11°.05 (considering that arcsin (i?ec/7?®) = 33°.11 and arcsin (i?jc/7?©) = 
11°.05). 

Thus, in this scenario, the DM particles move in each Earth’s layer 
with a mean free path \l, that mainly depends on the interaction cross 
section (see eq. 8). The number of interactions in each layer, Uhn, has 
been estimated as: 


- Case 1, high interaction cross sections, d© > 50A^: nhu is relatively 
high and follows a gaussian distribution with mean value and vari¬ 
ance equal to di/Xi] 

- Case 2, small interaction cross sections, d® < bOA®: a step-by-step 

approach has been adopted in the simulation. The path between two 
consecutive interactions, Xk, follows the distribution it 

^This equation holds in the limit where the form factor can be neglected. For high velocity, high 
mass DM candidate crossing the Inner Core it is only an approximation. For simplicity we do not 
consider further this issue here, while the form factor is considered in obtaining the energy loss (see 
later). 



can be used to propagate the particle within the layer as long as 
Xk < dL. 

In the considered scenario the DM candidate particles interact via SI 
elastic scattering on nuclei; thus, their energy-loss for each interaction is 
given by the induced nuclear recoil energy: 



( 10 ) 


where Em = mDMv‘^/‘^ and v are the DM energy and the velocity before 
the interaction, 6* is the angle of diffusion in the center of mass and r is 
a kinematic factor: 

AmDMrrii 

^ = 7-T-7^- 

[moM + rriiY 

The interaction cross-section is given by [18]: 

d(TDM,if \ d(j£iM,i f n\Tp2/rp \ 

Er) = {v, 0)E- [Er), (12) 

with E‘1{Er) nuclear form factor. In addition, for a given velocity v. 

do'DM,i, daDM,i dVt 

^' ie;' 

with (also see eq. 10): 

d(JDM,i _ 0'DM,i dER _ rEin , . 

dn ~ Att ’ YeV ~ 1^' ^ ^ 

Hence, assuming the cross section scaling law given in eq. 7, the interac¬ 
tion cross-section can be written as: 

do'DM,i , \ CTDM,! j-,2 / 7-1 \ A 2 -C'2( rp \ 


Thus, after one interaction, such a DM particle loses an energy in the 
range 0 < Er < rEm with distribution given in eq. 15. In the following 
we assume for the nuclear form factor the Helm formula^ [53, 54]: 


3lt is important to remark that the spin independent form factor depends on the target nucleus 
and there is not an universal formulation for it. Many profiles are available in literature and whatever 
profile needs some parameters whose value are also affected by some uncertainties. The form factor 
profiles can differ - in some intervals of the transferred momentum - by orders of magnitude and the 
chosen profile strongly affects whatever model dependent results [25]. 
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where Tq = — 5s^, r* = 1.2A\^^ fm, s ~ 1 fm and = 2miEji. 

In the simnlation, once nhu has been evalnated, the energy loss of the 
considered DM particle, is estimated for each interaction and the 
ontpnt velocity is: 





2-E r 
tudm 


(17) 


Therefore, the net effect is a modihcation of the velocity distribntion, 
fiah{y,t\mDM,o'n)] Fig. 3 sliows some examples of the obtained velocity 
distribntion for a detector located at LNGS. 
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Figure 3: Some examples of simulated fiab{'v,t\'mDM^<^ 1 %) in arbitrary units (a.u.) for 
some rriDM and tT„ values when taking into account the Earth Shadow Effect. In these 
plots: i) the considered Galactic Halo model is an isothermal sphere with vq = 220 
km/s; ii) the velocity distribution at GMST hour 8:00 (continuous black line) and 
20:00 (dashed - red on-line - line) - corresponding to the maximum and the minimum 
of the Earth Shadow Effect in case of a target-detector placed at LNGS. 
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4 The expected interaction rate 


In the SI coupling scenario, considered here, the DM candidates scatter 
off the nuclei in the detector. Their expected interaction rate as a function 
of the nuclear recoil energy, for a mono-atomic {i nucleus) detector 
is: 


dNj 

dEji 


Nj 


' '^min (Er) rriDM 


pv s.d(JDM,i 

Jlab{v] 


dE 


R 


V, 0)FpEg)dv, 


( 18 ) 


where Nt is the number density of the target nuclei in the detector and 
p = ^Po is the DM particles density in the galactic halo (.^ is the relative 
abundance and po the overall DM density in the galactic halo). The 
integral is calculated over all the possible DM particle velocities in the 
laboratory frame considering the distribution fiabiv) = fiab{v, t\mDM, o'n)- 


The minimal velocity providing E^ recoil energy is Vmin 



ER-mj 


The galactic escape velocity is included in the f{v) dehnition. Using eq. 
15, the expected rate can be rewritten as: 


dN. 


dE 


R 


= 


NtPo r.2 


-A: 


E/(E 


R) 


flab{v,t\mDM,0-n) 


rriDM 2 /i, 
dN' 


'^min (Er) 


dv 


dE 


( 19 ) 


R 


Generalizing to detectors with more than one kind of target nuclei (as e.g. 
in the case of the Nal(Tl) considered here), the expected experimental rate 
is: 


dEd^i 

where: 


dR dR' , I , 


'dE, 


det 


( 20 ) 


dR' 




dEfi^i 

[ G{E^„E' 


( 21 ) 


dN' 


Ki{E'\ER) {Er, t\mDM, CEn)dER 


dE': 


the G{Edet, E') kernel takes into account the detector’s energy resolution 
(generally through a gaussian convolution) and the Ki{E'\ER) kernel takes 
into account the energy transformation of the nuclear recoil energy in 
keV electron equivalent. For example the latter kernel can be written 
in the simplest case of a constant quenching factor, g*, as: Ki{E'\ER) = 
5{E' — qiER). For a discussion about the quenching factors see Refs. [5, 8]. 
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The expected differential rate - as well as fiab{v, t\mDM-, o'n) - depends 
on the time through three different effects: i) the time dependence of the 
Earth’s orbital motion velocity, Vrevit) (see eq. 1); ii) the time dependence 
of the Earth’s rotation velocity around its axis, Vrotit) (see eq. 1); iii) the 
possible Earth Shadow Effect which depends on Details about the 
Erst two effects, responsible of the model-independent annual and diurnal 
modulation of the DM signal rate, respectively, are reported in Ref. [12]. 
Following the same approach the expected rate in an energy interval 
can be written as^: 


Skit) = f 




dR' 


iEdeti t\rfl£)M 1 (^n)dE(i£^ 


Y dEdet / 

= ickn S'Q^^imDM.Ckn) + S'^f,i^DM,CFn)cOf^iuit-to)) ( 22 ) 

T ^d,ki^^DM 1 ^n) COs(caj.o^(t ^d)) T ^d,sh,ki^^DM i ^ni ^)j i 


where: i) Sq i^imDM,o'n) is the time independent component of the ex¬ 
pected signal; ii) Sl^i^imDM,o'n) is the annual modulation amplitude, 
u = 2 tt/T with T = 1 yr, to — June 2”^^; iii) S'^ j^irnDM, o'n) is the di¬ 
urnal modulation amplitude, Urot = ‘^^i/Td with Td = I sidereal day, td 
for the case of a detector at the Gran Sasso longitude ranges from 13.94 h 
to 14.07 h depending on the Vq value (see Ref. [12]); iv) S'^ ki^DM, o'n, i) 
(whose average value over Td is null) takes into account the signal vari¬ 
ation as a function of the sidereal time due to a possible Earth Shadow 
Effect. 

The ratio Rdy = S'^^^iniDM,(On)/S'n^^kiniDM,(yn) is model indepen¬ 
dent and it is Rdy — 0.016 at LNGS latitude; thus, considering the 
DAMA/LIBRA-phasel experimental result on the DM annual modula¬ 
tion, the expected fanS'dki'oo^DM,o^n) is order of 10“*^ counts per sidereal 
day per kg per keV (cpd^i^/kg/keV, hereafter) [12]. The reached ex¬ 
perimental sensitivity of DAMA/LIBRA-phasel [12] is not yet enough 
to observe such a diurnal modulation amplitude; in fact, in the (2-4) 
keV energy interval considered here, the experimental diurnal modula¬ 
tion amplitude from DAMA/LIBRA-phasel data is (2.0 ± 2.1) x 10“^ 
cpdsjrf/kg/keV [12] (< 5.5 x 10“^ cpd^j^/kg/keV, 90% G.L.). Thus, in the 
following we do not further approach it. 

Few examples of the S'^ j^it) behavior for tudm = 30 GeV and for dif¬ 
ferent values of the cross-section in the given framework are shown in Fig. 
4. For clarity, in the explicative case of (T„ = 10 pb and .^ = 1.1 x 10“® so 

“^Here only the first order terms are shown (i.e. the interference terms are omitted). 
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X 1000 



GMST (h) 

Figure 4: Examples of expected S'^ obtained following the approach 

described in the text. The Q/ quenching factor and the set A of the parameters’ values 
in presence of the channeling effect have been considered (see text); moreover, vq = 220 
km/s. The energy interval, considered here, is (2-4) keV and the mnu = 30 GeV. 
The obtained 5”^ Cn, t) are shown for the cases: i) an = 10 pb (continuous 

black line); ii) cr„ = 1 pb (dashed - red on-line - line); iii) cr„ = 0.1 pb (dotted - blu 
on-line - line); iv) cr„ = 0.01 pb (dot-dashed - green on-line - line). 


that ^an = 1.1 X 10“^ pb is compatible with the DAMA/LIBRA-phasel 
DM annual modulation result, the obtained amplitude of ^it) 

is of order of 3 x 10“^ cpd^j^/kg/keV. Such value can be studied in 
DAMA/LIBRA-phasel (see later). 


5 Data analysis 

The results, obtained by analysing in the framework of the Earth Shadow 
Effect the DAMA/LIBRA-phasel (total exposure 1.04 tonxyr) data, es¬ 
sentially depend on the most sensitive (2-4) keV interval; thus, this is the 
energy region considered here. 

In the present analysis, as in Refs. [5, 32], three possibilities for the Na 
and I quenching factors have been considered: Q/) the quenching factors 
of Na and I “constants” with respect to the recoil energy Er: q^a — 0.3 
and qi ~ 0.09 as measured by DAMA with neutron source integrated over 
the 6.5 — 97keV and the 22 — 330 keV recoil energy range, respectively 
[15]; Qii) the quenching factors evaluated as in Ref. [55] varying as a 
function of Er, Q///) the quenching factors with the same behaviour of 
Ref. [55], but normalized in order to have their mean values consistent 
with Q/ in the energy range considered there. 

Another important effect is the channeling of low energy ions along 


13 







axes and planes of the Nal(Tl) DAMA crystals. This effect can lead to 
an important deviation, in addition to the other uncertainties discussed 
above. In fact, the channeling effect in crystals implies that a fraction 
of nuclear recoils are channeled and experience much larger quenching 
factors than those derived from neutron calibration (see [5, 30] for a dis¬ 
cussion of these aspects). Since the channeling effect cannot be generally 
pointed out with neutron measurements as already discussed in details 
in Ref. [30], only modeling has been produced up to now. In particular, 
the modeling of the channeling effect described by DAMA in Ref. [30] 
is able to reproduce the recoil spectrum measured at neutron beam by 
some other groups (see Ref. [30] for details). For completeness, we men¬ 
tion an alternative channeling model, as that of Ref. [56], where larger 
probabilities of the planar channeling are expected. Moreover, we men¬ 
tion the analytic calculation claiming that the channeling effect holds for 
recoils coming from outside a crystal and not from recoils produced inside 
it, due to the blocking effect [57]. Nevertheless, although some amount 
of blocking effect could be present, the precise description of the crystal 
lattice with dopant and trace contaminants is quite difficult and analyt¬ 
ical calculations require some simplihcations which can affect the result. 
Because of the difficulties of experimental measurements and of theoret¬ 
ical estimate of this channeling effect, in the following it will be either 
included or not in order to give idea on the related uncertainty. 

Thus, the data analysis has been repeated in some discrete cases which 
allow us to account for the uncertainties on the quenching factors and on 
the parameters used in the nuclear form factors. The hrst case (set A) is 
obtained considering the mean values of the parameters of the used nu¬ 
clear form factors (see above and Ref. [25]) and of the quenching factors. 
The set B adopts the same procedure as in Refs. [20, 21], by varying (i) 
the mean values of the measured ^^Na and ^^^I quenching factors up to 
+2 times the errors; (ii) the nuclear radius, r^, and the nuclear surface 
thickness parameter, s, in the form factor from their central values down 
to -20%. In the last case (set C) the Iodine nucleus parameters are hxed 
at the values of case B, while for the Sodium nucleus one considers: (i) 
^^Na quenching factor at the lowest value measured in literature; (ii) the 
nuclear radius, r*, and the nuclear surface thickness parameter, s, in the 
SI form factor from their central values up to -|-20%. Finally, three values 
of Uq have been considered: i) the mean value: 220 km/s, and ii) two 
extreme cases: 170 and 270 km/s. 

Because of the large number of the needed simulations, the mass of the 
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DM candidate and of the cross section on nncleon have been discretized 


as in the following: six niDM (5 GeV, 10 GeV, 30 GeV, 60 GeV, 100 GeV 
and 150 GeV) and eight cr„, (10 pb, 5 pb, 1 pb, 0.5 pb, 0.1 pb, 0.05 pb, 
0.01 pb and 0.005 pb). 


% 0.04 



-0.04 


ill 


0 2 4 6 8 10 12 14 16 18 20 22 24 

GMST (h) 


Figure 5: Experimental model-independent diurnal residual rate of the single-hit scin¬ 
tillation events measured by DAM A/LIB RA-phasel in the (2-4) keV energy inter¬ 
val (crosses) as function of the sidereal time [12], superimposed to two examples of 
expectations obtained by MC simulation. The template curves are obtained in the 
same scenario as in Fig. 4, considering = 1.1 x 10“^ pb (compatible with the 
DAMA/LIBRA-phasel DM annual modulation result): i) the dashed (red on-line) 
line is obtained for cr„ = 10 pb (and, thus, ^ = 1.1 x 10“®); ii) the dotted (blu on-line) 
line is obtained for cr„ = 0.1 pb (and, thus, ^ = 1.1 x 10“®). The latter is compatible 
with the absence of diurnal rate variation in DAMA/LIBRA-phasel, while the former 
is not. 

The expectations are compared with the experimental model- 
independent dinrnal residnal rate of the single-hit scintillation events, 
measnred by DAMA/LIBRA-phasel in the (2-4) keV energy interval, 
as fnnction of the sidereal time (see in Fig. 2 of Ref. [12]). Two examples 
of expected signals are reported in Fig. 5. Here, the nsed sidereal time bin 
is 1 honr (24 time bins in the sidereal day) and the experimental residnals 
are: ± ad,i- We compnte the quantity: 



( 23 ) 


xl - 2Bi + Ae, 


where 
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The S'^ AniDM, CTnyU) have been evaluated for each set of parameters 
described above. 

The of OQ- 23 is function of only one parameter, Since the data 
do not show the presence of signihcant diurnal variation in the counting 
rate as already described in Ref. [12], only upper limits for ^ are allowed, 
once given rriDM and an- We can dehne: 


Ax^{e} = - xo- 


The Ax^ is a with one degree of freedom and is used to determine the 
upper limit of ^ parameter at 2a C.L.. 

-2 -2 




(pb) (pb) 


Figure 6: Examples of comparisons at 2a C.L. between allowed regions from the 
DM annual modulation results of DAMA/LIBRA-phasel (black continuous lines) and 
exclusion limits from the Earth Shadow Effect (dotted - red on-line - line). In both 
cases the quenching factors Q/, including the channeling effect, and the set B of the 
parameters have been considered (see text) for vg = 220 km/s and for two DM masses. 
(Left): muM = 10 GeV, the upper limits on ^ do not constrain the results of annual 
modulation. {Right)'. rnoM = 60 GeV, the upper limits on ^ do exclude the band 
with an > 0.05 pb and ^ > 10“^ for the considered model framework. The combined 
allowed regions are reported as shaded - green on-line - area. 


Two examples to describe the followed procedure are reported in Fig. 
6 , where the excluded regions (above the dotted lines) in the f vs cr„ plane 
for the cases of tudm = 10 and 60 GeV are shown as obtained on the basis 
of the Earth Shadow Effect in the given model framework. 
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The upper limits on ^ can be compared with the positive results from 
the DM annual modulation signature achieved by DAMA^. In particu¬ 
lar, DAMA/LIBRA-phasel reports an annual modulation amplitude in 
the (2-4) keV energy interval: = (0.0167 ± 0.0022) cpd/kg/keV, 

corresponding to 7.6cr C.L. [4]. 



Figure 7: Examples of the mean values of the allowed region of ^ as function of 
cr„ and rriDM, represented as an allowed surface (see text). The plots have been 
obtained for vq = 170 km/s in the following scenarios: a) the quenching factors Q/, 
without channeling and marginalizing over the parameters sets A, B, C; b) as in case 
a) including the channeling effect; c) the quenching factors Q//; d) the quenching 
factors Q///. We note that the “thickness” of the allowed regions around the surfaces 
is < ±30%; therefore, for simplicity it is not represented in these Figures. Finally, we 
recall that other uncertainties not considered in the present paper are present and can 
extend the result. 

Here for each set of parameters described above, one can evaluate (see 
e.g. eq. 22) the allowed values as: 

oexp 

__ 

^rn,{2-4)kev(^DM, ^n) 

^We recall that DAMA/LIBRA and the former DAMA/Nal have cumulatively reached a model 
independent evidence at 9.3<t C.L. for the presence of DM particles in the galactic halo on the basis 
of the exploited DM annual modulation signature [4]. 



( 24 ) 
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(a) 


(b) 






Figure 8: Examples of the mean values of the allowed region of ^ as function of (t„ 
and muM, represented as an allowed surface (see text). The plots have been obtained 
for vq = 220 km/s in the same sets of parameters as in Fig. 7. See text. 


(a) 


(b) 






Figure 9: Examples of the mean values of the allowed region of ^ as function of (t„ 
and moM, represented as an allowed surface (see text). The plots have been obtained 
for Vo = 270 km/s in the same sets of parameters as in Fig. 7. See text. 
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This corresponds, once including the experimental uncertainties on 
to a band in the ^ vs an plane (within the continuous solid line). In 
Fig. 6 such bands at 2a C.L. are reported. One can see that for the 
scenario considered there and for ttidm = 10 GeV the upper limits on 
^ do not constrain the results of the DM annual modulation. On the 
contrary, for uidm = 60 GeV the upper limits on ^ do exclude the band 
with an > 0.05 pb and ^ > 10“^. The shaded bands in Fig. 6 corresponds 
to the allowed regions in the ^ vs plane, for the given moM, from the 
combined analyses of the DM annual modulation result and of the Earth 
Shadow Effect in the considered framework. 

Finally, for each considered set of parameters the three-dimensional 
allowed region - calculated as described above - in the parameter’s space: 

an, moM, is depicted as a surface in Fig. 7, 8 and 9, for vq equal to 170, 
220, 270 km/s, respectively. We note that the “thickness” of the allowed 
regions around the shown surfaces is < ±30%; therefore, for simplicity it 
is not represented in these Figures. 

Finally, we recall that other uncertainties not considered in the present 
paper are present. For example, including other possible halo models 
sizeable differences are expected in the results as shown e.g. in Refs. 
[5, 22]. 


6 Conclusions 

The Earth Shadow Effect has been investigated in a given framework 
considering the model independent results on possible diurnal varia¬ 
tion of the low-energy rate of the single-hit scintillation events in the 
DAMA/LIBRA-phasel data (exposure: 1.04 tonxyr) reported in Ref. 
[12]. For the considered DM candidates having high interaction cross- 
sections and very small halo fraction the obtained results constrain at 
2a G.L., in the considered scenario, the an and niDM parameters (see 
Figs. 7, 8 and 9) when including the positive results from the DM annual 
modulation analysis of the DAMA/LIBRA-phasel data [4]. For example, 
in the considered scenario for quenching factors Q/ with channeling ef¬ 
fect, B parameters set, Uq = 220 km/s and uidm = GeV, the obtained 
upper limits on f do exclude an > 0.05 pb and f > 10“^. When also 
including other uncertainties as other halo models etc. the results would 
be extended. 
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